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ABSTRACT 

We have observed seven ultracompact H II regions in hydrogen recombination 
lines in the millimeter band. Toward four of these regions, there is a high velocity 
(full width to half maximum 60-80 km s^^) component in the line profiles. 
The high velocity gas accounts for 35%-70% of the emission measure within 
the beam. We compare these objects to an additional seven similar sources we 
have found in the literature. The broad recombination line objects (BRLOs) 
make up about 30% of all sources in complexes containing ultracompact H II 
regions. Comparison of spectral line and continuum data implies that the 
BRLOs coincide with sources with rising spectral indices, >0.4 up to 100 GHz. 
Both the number of BRLOs and their frequency of occurrence within H II region 
complexes, when coupled with their small size and large internal motions, mean 
that the apparent contradiction between the dynamical and the population 
lifetimes for BRLOs is even more severe than for ultracompact H II regions as 
a whole. We evaluate a number of possible models for the origin of the broad 
recombination line emission. The lifetime, morphology, and rising spectral 
index of the sources argue for photoevaporated disks as the cause for BRLOs. 
Existing models for such regions, however, do not account for the large amounts 
of ionized gas observed at supersonic velocities. 

Subject headings: ISM:H II Regions ISM:Jets and Outflows Radio Lines:Stars 
Stars:Mass Loss 
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1. Introduction 

The formation, behavior, and evolution of O stars influence many aspects of the 
interstellar medium in the Galaxy. O stars play a central role in the evolution of galactic 
abundances, influence the structure of the interstellar medium through supernova explosions, 
ionization, and mass-loss, and determine the appearance of the spiral arms in the visible 
and the radio through both thermal and non-thermal emission. The pre-main sequence 
evolution of luminous O stars takes place while they are still shrouded within their natal 
clouds. Therefore, we usually study the surrounding H II regions and dust clouds, rather 
than the stars themselves, and carry out these studies at far-IR and radio wavelengths. 

Much of the gas in the H II regions that surround newly formed O stars is involved 
in supersonic motions of the order of 20 km s~^, usually ascribed to "turbulence". The 
smallest and presumably youngest H II regions (ultracompact, or UC H II regions) typically 
have linewidths of 25-35 km s~^ in high frequency centimeter wave recombination lines 
(RRL's, Garay 1989, Churchwell et al. 1989, Afiierbach et al. 1996). There is considerable 
observational evidence for even higher velocity motions of various sorts. Massive O stars 
have winds with terminal velocities of 1-2x10^ km s~^ and mass loss rates of up to a few 
10~^ Mq per year (Chlebowski & Garmany 1991). Many high luminosity and presumably 
young stars embedded deep within molecular clouds show evidence for ionized circumstellar 
envelopes or for winds with terminal velocities ~200 km s~^ (Simon et al. 1981, 1983, 
Persson et al. 1984, Garden & Geballe 1986, Bunn, Hoare, & Drew 1995). In some more 
evolved and optically accessible H II regions, there are small amounts of gas moving at 
100-200 km s~^ with respect to the bulk velocity of the H II region (Meaburn & Walsh 
1981, 1986). 

There is increasing evidence that a (substantial) minority of ultracompact H II regions 
emit broad radio recombination lines. (In this paper, we refer to these sources as broad 
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recombination line objects, BRLOs). The presence of this RRL component imphes that the 
regions contain gas moving at velocities up to ±60 km s~^ from the rest velocity of the cloud 
(Altenhoff, Strittmatter, & Wendker 1981, Zijlstra et al. 1990, Gaume, Fey, & Claussen 
1994, De Free et al. 1994, 1996). In the UC H II regions, unlike the more evolved H II 
regions, the high velocity gas may account for a significant fraction of the total emission 
measure. If high velocity emission arises in the very compact parts of the UC H II region 
complexes where the emission measures are high, the gas may only become detectable in 
observations at the shortest centimeter wavelengths where these compact sources begin 
to become optically thin and where non-LTE effects become negligible. Observations at 
shorter wavelengths could therefore provide a better way to search for the broad features. 
We have investigated the incidence and nature of the broad recombination line features by 
surveying a sample of bright UC H II regions in 3, 2, and 1.3 mm recombination lines. This 
paper presents results and analysis of this survey and a discussion of the origin of the broad 
line emission. 

2. Observations 

We made all the observations we present here with the IRAM 30 meter telescope at 
Fico de Veleta (Spain). The first set of observations took place in 1992 February and the 
second set in 1993 March. On the fist observing run, we observed three recombination 
lines simultaneously at wavelengths of 3, 2 and 1.3 mm, and on the second run, we 
observed several combinations of two such lines (see Table |^). The telescope beam at 
these wavelengths was, respectively, 24, 17 and 11 arc seconds. Table 1 lists the sources, 
source positions, dates of observations, and observed transitions. The receivers were SIS 
mixers, tuned to single sideband (image rejections between 7 and 10 dB), with typical 
system noise temperatures of 250, 350 and 650 K at 3, 2 and 1.3 mm. We used two 512 
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channel filterbanks with 1 MHz resolution and a 430 channel Acousto-Optical Spectrometer 
with a resolution of 1.1 MHz as spectrometric backends. In order to insure good baseline 
subtraction, we nutated the telescope subreflector at 0.25 Hz between the source and a 
reference position 240 arc seconds away. While some of the sources have extended radio 
continuum emission, published maps at centimeter wavelengths show that an angular 
displacement of this size should result in contamination of the reference beam by narrow 
recombination line emission at only a very low level. Under good weather conditions, 
use of the nutating subreflector allows one to measure the line and continuum emission 
simultaneously. The weather conditions during the 1993 observations were good enough 
to permit us to measure the continuum emission at 3 mm. We calibrated the data by 
observing the atmosphere and cold and hot loads at known temperatures. All the intensities 
are in units of antenna temperature corrected for forward scattering and spillover as well as 
for atmospheric transmission, T^, and the conversion factor from antenna temperature to 
flux density at 3 mm is 6 Jy/K. 

Almost all of the ultracompact H II regions studied as part of the present work lie 
close to or within dense, high column density molecular cores. These cores produce very 
rich millimeter wave line spectra. Confusion from these molecular line emitters within 
the beam is an important limiting factor when deriving quantitative results from the 
millimeter recombination line data. The presence of the molecular lines does not, however, 
impede our ability to detect high velocity recombination line emission. The molecular lines 
are, in general, quite narrow (5-8 km s~^) and the broad components of adjacent An=l 
recombination transitions with similar n should be similar in strength and shape but overlap 
with completely different sets of molecular features. 
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3. Results and Analysis 

Table |^ presents the best fits of one or two Gaussians to profiles of selected lines 
in each source in our survey. In all cases, the lines selected from among the 1.3-3 mm 
Hna transitions observed (Table jT]) are the 3 mm transitions, H39a and H41a. Since our 
observed positions coincide with or lie close to hot, dense, high column density molecular 
cores (cf. Plume, Jaffe, & Evans 1992, Plume et al. 1997), the high density of molecular 
features at 2 mm and 1.3 mm makes the hydrogen recombination lines in these bands more 
prone to confusion. The 3 mm continuum flux densities measured toward the positions 
given in Table |l| also appear in Table ^ Table |^ gives the electron temperatures, derived 
from the ratio of the total integrated H41a line intensity (including both line components 
if more than one is present) to the measured continuum adjacent to this 3 mm line, 
assuming LTE conditions. It is interesting to note that failure to notice the presence of a 
relatively weak broad component in recombination line spectra, a component which can be 
missed in observations with moderate sensitivity, can lead one to overestimate the electron 
temperatures derived from the line-to-continuum ratios. The LTE electron temperatures 
we derive from our data all lie within ±30 % of the canonical 10,000 K temperature of H II 
regions. 

3.1. Broad Recombination Line Objects 

The millimeter recombination line shapes for four of the seven sources in our sample 
(NGC 7538, W49N, G34.3+0.2, and S106) deviate substantially from the simple 25-35 km 
s~^ wide Gaussian profiles that are typical of ultracompact H II regions. Fits of single 
Gaussians to these profiles left systematic residuals at |y| >20 km s~^ from the line center. 
In all cases, we significantly reduce the residuals when we fit the observed spectra with 
two Gaussian profiles, one with the normal UC H II region recombination line width and 
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another broad component which is 2-3 times wider. We present here our analysis of these 
four sources. 

Toward most of the sources in our sample, the 30 m beam contains very complex 
radio continuum structure that contributes to both the normal and the broad spectral 
components. So far, there are no maps of recombination lines at millimeter wavelengths 
with sufficient angular resolution and sensitivity to identify the sources responsible for 
the broad features. To identify the sources giving rise to these features, we combine our 
millimeter recombination line data with existing high angular resolution maps of continuum 
and recombination line emission, mainly observed at centimeter wavelengths. To aid in 
comparisons with these maps, we estimate the amount of continuum emission "required" 
to produce each line component by assuming that the observed 3 mm line and continuum 
emission is optically thin and that the millimeter recombination lines arise from gas where 
the electron temperature equals the temperature derived from the total line and continuum 
emission (under the assumption of LTE conditions) from each source (Table ^. Table § 
gives our estimate of the continuum associated with the broad and the normal components 
based on these assumptions. The gas that produces the broad line component accounts for 
35%- 70% of the 3 mm continuum flux and therefore at least that fraction of the emission 
measure in these four objects. 

3.1.1. NGC 7538 

We have observed six millimeter wave hydrogen recombination lines toward NGC 7538. 
Our beam center lies at the position of the infrared source and double radio continuum 
source NGC 7538 IRSl (Campbell 1984). A free fit of a single Gaussian profile to the H41a 
line underestimated the brightness both at the peak of the line and at high velocities. The 
left panel of Figure |1] shows the results of a fit which we forced to match the area around 



the peak by minimizing the least squares deviation between a single Gaussian profile and 
the observed line over a range of ±15 km s~^ about the peak of the line. The lower part of 
the left panel shows that, while this profile accounts well for the inner part of the observed 
line, there are substantial residuals to this fit, in particular at higher velocity. The right 
panels of Figure |l] show a fit of two Gaussians to the same profile and the residuals resulting 
from this fit. The line parameters we derive in an independent fit of two Gaussians to the 
H39a line agree well with the H41a results (Table Q). For both hydrogen lines, the normal 
line has a width of around 25 km s~^ while the broad component has a width of 70 km s^^. 

High resolution 1.3 cm recombination line observations toward NGC 7538 IRSl (Gaume 
et al. 1995) reveal extremely broad emission (full width to zero power of around 200 km s~^ 
in the northern and southern hypercompact lobes and an absence of line emission in between 
the lobes. The discrepancies between the linewidths of the centimeter and millimeter 
wavelength recombination lines are similar to those found in MWC 349. As in this latter 
source, the broad H66a lines in NGC 7538 IRSl, which show noticeable non-Gaussian 
profiles with fiat-topped shapes, are very likely affected by non-LTE excitation (see 
Martin-Pintado et al 1993). 

For NGC 7538, our beam contains the sources IRSl, IRS2 and IRS3 and the total 
continuum flux density at 3 mm is in good agreement with the fluxes derived by Akabane 
et al. (1992) from interferometric measurements. The 3" resolution 3 mm continuum map 
(Akabane et al. 1992) shows that IRS3 makes a negligible contribution to the millimeter 
wave continuum flux, IRS2 contributes 1 Jy, and IRSl contributes 1.6 Jy. The values for 
the continuum fluxes for the two latter sources are consistent with an attribution of the 
normal spectral component to IRS2 and the broad component to IRSl, the same sense in 
which the origin of the normal and very broad 6 cm recombination line emission is directly 
attributable to the respective sources (Gaume et al. 1995). 
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3.1.2. W49N 

We observed six millimeter wave hydrogen recombination lines toward this source 
(Table ||). Our beam center lies toward the H2O maser center, about 1" north of the 
ultracompact 15 GHz continuum source Gl (Dreher et al. 1984). The upper left panel 
of Figure § shows a single Gaussian fit to the H41q: line made following the procedure 
described for the single Gaussian fit to the NGC 7538 profile. The residuals to this fit (lower 
left) show a systematic excess in the wings of the line. The upper right panel of Figure 
shows a two Gaussian fit to the same line. The minimization results in a component 
with a velocity width of 31 km s~^ and a second component 71 km s~^ wide. Over the 
velocity interval occupied by the hydrogen line, the residuals to the two Gaussian fit at the 
central position are flat. We also fit two Gaussians to the He41a; line at this position. We 
allowed the temperatures of the two components to vary but fixed the widths of the two He 
components at the values of the H41a widths and the velocity offset from the corresponding 
H41a components at -122.2 km s^^. This fit (Figure |^ (lower right)) shows that a normal 
plus a broad component can account for the shape of the helium line as well. The strengths 
of the normal and broad He components are respectively 9% and 8% of the strengths of the 
corresponding hydrogen emission components. A five point map of the H41a transition at 
12" intervals is consistent with an unresolved source as the origin of both the normal and 
the broad component. 

The upper panel of Figure |^ superposes spectra of H39a and H41a toward W49N. 
The wings of the two recombination lines match very closely, but the shapes of the inner 
portions of the lines differ considerably. A triple Gaussian fit to the H39a spectrum fixing 
the velocities and widths of the two broader components to match the values derived from 
the H41a spectrum shows that there is a third feature with a full width to half maximum 
of 13.3 km s~^. This narrower feature has a width similar to the width of molecular lines 
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in this source (Serabyn, Guesten & Schulz 1993). The third hne has an observed velocity 
of -7.5 km s~^ with respect to the H39a frequency. The only molecular transition with a 
rest frequency close to that of the H39a line (106737.365 MHz) is the 23 — I2 transition of 
^^SO (i^o=106743.365 MHz). For this line, the observed velocity corresponds to 'Vlsr= 6.6 
km s~^ in its own frame. This velocity is close to the typical molecular velocity of the W49 
core. The lower panel of Figure ^ shows an overlay of the observed H39q; and H41a profiles 
where we have added the fitted ^''SO feature to the H41a data to illustrate the similarity of 
the shapes of the recombination line profiles themselves. 

There is some evidence for high velocity recombination line emission in previous 
millimeter observations of W49N. The wings seen in Figure ^ are also apparent to the 
trained eye in the H40q; spectrum in Figure 3 of Gordon (1989). High angular resolution 
observations of the H66a and H52a transitions at 1.3 cm and 0.7 cm, respectively (De Free, 
Mehringer, & Goss 1997), show that sources A and B, both of which lie within our beam, 
produce lines with a full width to half maximum of 50-60 km s~^. Since the line fills the 
spectrometer band for these sources, however, the observations do not exclude the presence 
of higher velocity emission consistent with our millimeter measurements. 

Determining the source or sources in W49N responsible for the broad line emission is 
complex because of the large number of continuum sources in our beam. The H41a beam 
contains the sources A-G2. We can use high angular resolution maps of the RRL and 
continuum emission (De Free et al. 1997) to establish which source may be responsible for 
the broad component. As previously mentioned, sources A and B have the broadest H66q; 
and H52a profiles. These two sources have power law radio continuum spectra with spectral 
indexes of 0.6 and 0.9, respectively. Extrapolating to 3 mm, A and B would each contribute 
around 2 Jy. Given the 3 Jy required to provide the source of the broad line emission (see 
Table ^ , these two stellar wind sources are viable candidates for the origin of the broad 



-11- 



lines. However, source A is located at the edge of the beam and the map of the H39a hne 
indicates a point source located closer to B than to A. Like in NGC 7538 IRSl, the broad 
component seems to be associated with compact sources with rising continuum spectra. 

3.1.3. G34.3+0.2 

We observed six millimeter wave hydrogen recombination lines toward this source 
(Table |l|). The center position for these observations lies about 3" south of the 15 GHz 
continuum peak (Garay, Rodriguez, & van Gorkom 1986). The left panel in Figure ^ shows 
a single Gaussian fit to the H41a transition toward G34.3+0.2 (see the discussion of NGC 
7538 for a description of the fitting procedure). The residuals (lower left) show significant 
excess emission at high velocities, particularly on the blue side of the line. The right panel 
of Figure § shows a fit using two Gaussians. Superposed on the residuals of this fit is a 
fit to the helium line made by holding the relative velocity to the corresponding hydrogen 
component fixed and setting the widths of the two helium components equal to those of the 
hydrogen components. Figure |^ shows a superposition of the H41a and H39a lines toward 
G34.3+0.2. The figure also shows the sum of two Gaussians with parameters averaged 
between those from the H41a and H39a line fits. The compromise Gaussians have widths 
of 32 and 75 km s^^. The shape of the two hydrogen lines agree well, except at velocities 
where narrow molecular lines from the hot core contaminate the profiles. 

Toward G34.3+0.2, our 30 m beam at 3 mm contains three radio continuum sources. 
Components A and B are compact objects with spectra that rise with frequency at least 
to 90 GHz (Gaume et al. 1994). Component C is a classic cometary HIT region (Reid 
& Ho 1985; Garay et al. 1986). At centimeter wavelengths, where sources A and B are 
weak and opaque, all of the recombination line emission comes from source C. The 2 cm 
results (Garay et al. 1986) indicate a relatively uniform line width of around 50 km s~^ 
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across source C The more sensitive 3 cm data (Gaume et al. 1994) show hne widths which 
increase from 25-45 km s^^ in the north and west to about 80 km s~^ in the southeast of 
component C (Gaume et al. 1994). 

Our total continuum emission measured at 3 mm is ~7 Jy, in good agreement with the 
6.6 Jy we extrapolate from the 6 and 2 cm measurements of Benson & Johnston (1984) by 
assuming that source C is optically thin and ultracompact sources A and B are optically 
thick with spectral indexes of ~0.9 (Gaume at al. 1994). Based on this extrapolation, the 
compact sources A and B, despite their rising "stellar wind" spectra, can only account for 
0.8 Jy of the 7 Jy flux density at 3 mm. Both the broad and the narrow line component 
therefore should arise mainly from the cometary H II region. The radial velocities and 
widths of both the broad and narrow line components in the H41a spectrum lie within the 
range of one or another centimeter RRL measurement toward the parts of source C that 
emit most of the flux (Garay et al. 1986, Gaume et al. 1994). The results obtained by the 
centimeter RRL observers, however, are not consistent with each other. 

3.1.4. S106 

S106 is less luminous than the other BRLOs in our sample. We selected this 
source for investigation because of the numerous observations of broad near-IR hydrogen 
recombination lines toward S106 IR (Persson et al. 1984, Felli et al. 1985, Garden & 
Geballe 1986, Drew, Bunn, & Hoare 1993). This IR source lies at the center of an obscuring 
dust disk (Eiroa, Elsaesser, & Lahulla 1979) and at the edge of our beam, 13" northeast 
of the position listed in Table |l[ Coincident with the IR source is a compact radio source 
with a spectrum rising as v^'^ from 5 to 23 GHz (Snell & Bally 1986, Felli et al. 1985). 
Very high angular resolution 5 GHz observations show that this central source is elongated 
along the equator of the extended bipolar optical and radio nebula (Hoare et al. 1994, Felli 
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et al. 1984). There is a strong submillimeter dust continuum source close to the position 
we observed (Richer et al. 1993). The difference between the observed position and the 
IR source position and the presence of relatively strong dust emission make modeling this 
source problematic. 

For SI 06, the broad component cannot arise only from the compact continuum 
source with the rising spectrum since the extrapolation of the continuum spectrum of 
this source to 3 mm implies that it will only account for 42 mJy at this wavelength. 
Therefore, the high-velocity recombination line emission is associated with the bipolar 
nebula. Like G34.3+0.2, this source appears to be more evolved than the UC HII regions 
in NGC 7538 IRSl and W49N. Long slit spectroscopy in the optical shows that S106 has 
the velocity structure of a bipolar outflow (Solf & Carsenty 1982). There is Ha emission 
throughout the southern lobe of the extended bipolar source. Ha has a total extent of over 
100 km s~^ at the position we observed at 3 mm. Farther south in this lobe, the Ha lines 
split into two components with typical separations of 50 km s~^ (Solf & Carsenty 1982). 

3.2. Narrow Line Sources 

3.2.1. Orion/ IRc2 and Orion H2 Peak 1 

Figure ^ shows an H39a spectrum toward H2 Peak 1 {left) and an H41a spectrum 
toward IRc2 {right) in Orion. In both cases, the residuals below each spectrum demonstrate 
that a single Gaussian fits the line profile well. Hasegawa & Akabane (1984) have reported 
high velocity emission toward Peak 1 in the H51a line. The reported high velocity emission 
gave rise to the suggestion that such emission could come from gas ionized by the radiative 
precursor of a dissociative shock in the outflowing molecular gas (McKee & Hollenbach 
1987). Our H39a observation at the position of Peak 1 excludes the presence of a broad 
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component at a level about a factor of five lower than the 0.1 K implied by the putative 
H51a detection. Our H41a observations toward IRc2 also show no evidence for emission 
beyond that arising from a single 23 km s~^ wide Gaussian component. There is therefore 
no observational evidence for shock- ionized gas in the region around Orion/IRc2. It 
would be better to search for recombination lines from shock-ionized gas where little or 
no confusing photoionized gas is present. Possible targets might include the high velocity 
shocked regions in powerful outflows associated with less luminous sources like IRAS 03282 
and L1448 (Bachiller 1996). 

3.2.2. Gl 0.6- 0.4 

GlO.6-0.4 is a bright UC HII region situated in a molecular core. Observations of the 
radio continuum at 6, 2, and 1.3 cm show several emission knots within 5" (Ho & Haschick 
1981, 1986). We have observed two lines toward this source, 1135a; and H41a (Figure 0). 
The single Gaussian fit to the less confused II41a line (Table 0) undershoots the peak of 
the line and overshoots in the wings. This behavior mimics that of the residuals of single, 
unconstrained Gaussian fits to the lines in G34.3+0.2 and W49, indicating the possible 
presence of a broad component. The signal to noise ratio of the GlO.6-0.4 spectra, however, 
is insufficient to warrant a fit of more than one Gaussian. 



3.2.3. DR21 

We observed the H30a, II35a, H36a, H39a, and H41a transitions toward the compact 
H II region in DR21. Figure || shows the II39q; line. The residuals in this figure indicate that 
a single Gaussiam with a width of 32 km s~^ can account for all of the fiux in the hydrogen 
line. This source has strong molecular and neutral atomic emission at high velocities (Bally 
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& Lada 1983, Russell et al. 1992) but does not have a corresponding ionized gas component 
in the millimeter recombination lines. 



4. Other Broad Recombination Line Objects 

A search of the literature shows that there are at least seven additional ultracompact 
HII regions with high velocity (AV>55 km s~^) emission in recombination lines at n<80, 
where pressure broadening is less likely to be important. These sources are: the four sources 
comprising SgrB2M F, G5.89-0.39, MWC 349, and K3-50A (see Table 1). 

As a whole, SgrB2M source F has broad H66a and H52a lines (Mehringer et al. 1995, 
De Free et al. 1996). Integrating over the source, the H52a line has a FWHM of 59 km 
s~^. The H66a; observations with the VLA resolve source F into four subsources. The line 
widths for these more compact objects range from 68 km s~^ to 79 km s~^. Each of the 
subsources has a rising continuum spectrum at 1 cm. 

G5. 89-0. 39 contains an ultracompact source (Source A) with broad RRL emission. 
Zijlstra et al. (1990) detected a normal (AV= 22 km s^^) and a broad (AV= 64 km s~^) 
H76a component toward this source. Their continuum measurements imply a spectral 
index of v^'^ from 20 cm to 6 cm and z/°'* from 6 cm to 2 cm. AfHerbach et al. (1996) 
measured line widths of 59 km s~^ and 61 km s~^ at H42a and H66q;, respectively. Their 
data imply a spectral index of z/°'^^ from 1.3 cm to 3 mm. Additional continuum data in 
the literature give somewhat inconsistent results (Acord, Churchwell, & Wood 1997; Acord, 
Walmsley, & Churchwell 1998; Wood and Churchwell 1989), but all of the results imply a 
rising spectrum throughout the centimeter band and continuuing to at least 7 mm. Acord 
et al. (1998) have used 3.6 cm continuum observations over a 5 year span to measure what 
they argue is a proper motion of a shell due to dynamical expansion. 
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MWC 349 is an isolated HIT region in Cygnus. This object is the archetypal ionized 
stellar wind source. Its spectrum rises as ly^'^ from cm wavelengths up to 60 fim. (Harvey, 
Thronson, & Gatley 1979). This behavior is characteristic of an isothermal wind expanding 
at constant velocity (see e.g. Wright & Barlow 1975). The typical line widths measured 
for the centimeter wave recombination lines (AltenhofF et al. 1981) and optical [N II] 
lines (Hartmann, Jaffe, & Huchra 1980) is ~ 60 km s^^. VLA observations of the II66a 
line show that the lines have nearly flat topped profiles with a width of around 100 km 
s^^, substantially larger than those measured from the Gaussian profiles of the 3 mm 
recombination lines (Mart in- Pint ado et al 1989). Modeling of the continuum and the line 
intensities of the H66q; lines indicates that the discrepancies are likely to be due to non-LTE 
effects in the centimeter recombination lines (see Martin-Pintado et al 1993). MWC 349 
shows a bipolar morphology (White & Becker 1985), but unlike K3-50A (see below), it does 
not show a bipolar distribution in the recombination lines (Martin-Pintado et al 1993); 
indicating that MWC 349 must be oriented perpendicular to the line of sight. As an object 
for study, MWC 349 has the advantage over the other BRLOs detected so far that it is 
isolated and the ionizing star can be seen in the visible and the near-IR. In some sense, 
MWC 349 may play the role for massive stars that field TTauri stars play in studies of low 
mass stars. 

At 2-6 cm, the elongated continuum source K3-50A dominates the radio continuum 
emission from compact radio sources in the K3-50 complex. Combining the continuum data 
of Roelfsema, Goss, & Geballe (1988) and De Pree et al. (1994), we estimate a 5-15 GHz 
spectral index for this source of u^'^. Roelfsma et al. (1988) detected a broad (AV= 66 km 
s~^) component to the HllOa line but attributed its breadth to pressure broadening. The 
more recent H76q; results (De Pree et al. 1994) show not only that the large line widths are 
a dynamical effect but also that there is clear evidence in this source for an ionized bipolar 
outflow. Toward the center of the source, the 2 cm recombination line has a width of 76 km 



-17- 

s~^. Three arc seconds north and south of this position, the width is only 35 km s~^ but 
the velocity changes from -11 km s~^ north of the peak to -42 km s~^ to the south. 

5. Discussion 

5.1. Observational Constraints on Possible Models 

In this section, we summarize the observations and key parameters derived from the 
observations for which any viable model for objects that produce the broad recombination 
lines (BRLOs) must account. 

5.1.1. Lifetimes 

BRLOs are common in UC H II region complexes. The lifetimes of these sources are 
much longer than their crossing times and the broad line emission phase must represent a 
substantial fraction of the lifetime of UC H II sources. 

One argument for the long duration of the broad line emission phase as a part of overall 
H II region evolution is the presence of broad millimeter and high frequency centimeter 
recombination lines toward ultracompact sources like NGC 7538 IRSl and MWC 349 and 
also toward more extended and presumably more evolved sources like K3-50 and S106. It 
is possible that there is a connection between the dominant broad line emission at early 
evolutionary stages and the low-level emission at high velocities seen in very evolved HII 
regions (Meaburn & Walsh 1981, 1986). The available sample of H II regions observed 
with high enough sensitivity to detect the broad recombination lines is small (Table ^. 
As a consequence, our estimate of the fraction of compact H II regions with broad lines 
is somewhat uncertain and may suffer from selection biases. We follow two approaches in 
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estimating the frequency of occurrence of these objects: We first determine the fraction 
of sources in UC H II region complex that show the broad lines, within each well-studied 
complex containing at least one BRLO. In NGC 7538, one out of three sources is a BRLO. 
In G34.3+0.2 one to three out of three is likely to be a BRLO and in W49, we find that 
two out of eight sources are probably responsible for the broad recombination line emission. 
From these data, we estimate that roughly around 30% of the compact sources in these 
complexes are BRLOs. 

As a second approach, we use a limited sample of compact H II regions (Afiierbach et 
al. 1996). Since we are only interested in the statistical properties of the sample and since 
the signal to noise in the recombination line data from this survey is lower than for any of 
the observations of individual sources, we set a somewhat lower threshold for BRLOs and 
include objects with RRL widths >45 km s~^. Using this criterion, 6 of the 23 sources in 
the sample have lines broad enough to demonstrate the presence of one or more BRLOs 
which also implies a frequency of incidence of roughly 30% for BRLOs in their sample. 

Table |^ gives a summary of the main characteristics of all the well-studied BRLOs 
(i.e. not including the statistical sample from the AfHerbach et al. survey). Most of the 
sources are very compact, with typical crossing times of less than 2x 10^ years for expansion 
velocities of around 50 km s~^. Following the same line of argumentation used by Wood 
& Churchwell (1989) for all compact HII regions, we find that, for such short dynamical 
times, one would expect to detect only ~15 objects of this kind in the whole galaxy, close 
to the number of objects so far detected in just a few regions. The high detection rate for 
BRLOs means that these objects cannot be expanding at the velocities measured by the 
recombination lines. 
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5.1.2. Continuum Properties 

Most of the resolved BRLOs have a bipolar morphology and almost all BRLOs have 
rising continuum spectra up to at least 1.3 cm (Table |]). Even for the most compact 
objects, the rising spectrum continues well into the millimeter wavelengths. Unfortunately, 
the complexity of most of the regions and the presence of large column densities of dust 
makes it very difficult to measure the continuum spectral indexes at short millimeter and 
sub-millimeter wavelengths. For the best case, MWC 349, the v^'^ continuum spectrum is 
observed up to 60 ^m. (Harvey et al. 1979). If other BRLOs share this characteristic, large 
emission measures lie hidden in the optically thick regions seen at centimeter wavelengths 
in most of the BRLOs in Table 3. As determined from the high- resolution VLA maps, all 
the BRLOs have emission measures of at least 10^ pc cm~^ and the number of the Lyman 
continuum photons needed to explain the radio data varies from 10^^'® s~^ for S106 to 
more than 10^*^'^ s~^ for the other sources. For the sources with continuum spectra that 
rise into the millimeter wavelengths, these values are clearly lower limits and the emission 
measure and the number of Lyman continuum photons increase with frequency at least 
with the same power law as the continuum spectrum. This indicates that the ionization 
of all the BRLOs, except for S106, will require stars with spectral types earlier than 06. 
The large emission measure indicates that the BRLOs must be ionized by a star or stars 
at the geometrical centers of the HII regions, within the optically thick cores observed at 
centimeter wavelengths. The central location of the ionizing stars is already apparent in 
MWC 349 and S106 where the geometry is clear because the ionizing stars are observed in 
the visible and the IR. 

The rising continuum spectral index of an opaque ionized source is a product of a 
combination of the radial density and temperature distributions of the ionized material. 
This index is nearly independent of the axial and azimuthal geometry (Schmid-Burgk 1982). 
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The BRLOs mapped with the VLA at high angular resolution show that the continuum 
and recombination line intensities are consistent with minor variations of the electron 
temperatures within the sources (see e.g De Pree et al. 1994 for K3-50, Gaunie et al. 
1995 for NGC 7538 IRSl; Martin-Pintado et al. 1993), suggesting that the HII regions 
are close to isothermal. Furthermore, for the typical spectral indexes measured for the 
BRLOs, between 0.5 and 1.5, the effects of moderate power law temperature gradients have 
minor effects (~20%) on the continuum spectral indexes; being independent of temperature 
gradients for spectral indexes of ~0.6 (Cassinelli & Hartmann 1977). Therefore, we will 
consider that the spectral indexes given in Table 3 are basically related to the density 
structure within the BRLOs. The density structure varies from nocr~^ to nocr^^ for spectral 
indexes between ~0.6 and 1.3 respectively (Olnon 1975; Panagia & Felli 1975; Wright 
& Barlow 1975). For such steep density gradients, the HII region is far from pressure 
equilibrium and the ionized material must be flowing away from the dense region near the 
star. From mass conservation, the density structures derived from the continuum spectral 
indexes, imply that only in a small fraction of BRLOs (those with spectral indexes of ~ 0.6) 
is the ionized gas outflowing at constant velocity, and in most of the BRLOs the ionized 
gas is accelerating. In the most extreme case (spectral indexes of 1.3), the velocity must be 
proportional to the radius. 

5.1.3. Line Widths 

The measured linewidths of the BRLOs (Table ^) fall within a remarkably small range. 
Our best estimate for the completeness range for our data is AV(FWHM) =45-100 km s~^. 
Confusion with narrow line emission from the same complex imposes the lower bound. We 
estimate the upper bound from our observed spectra. The upper bound results from a 
combination of the diminuition of peak intensity with increased width and confusion with 
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the helium lines and with low-level molecular line emission. From our own observations, 
however, the widths of the broad recombination line components cluster strongly around 75 
km s~^. The results in Table ^ show that even sources observed at centimeter wavelengths 
with interferometers have broad recombination lines with widths very close to 75 km s~^. 
If further observations bear out this consistency in the velocity spread of the emission, the 
explanation of this result must form an important part of any theory of BRLOs. 

5.2. Possible Models 

In this section, we summarize the various possible models for UC H II regions in the 
context of the our new data and other recently published data on broad recombination 
line emission. We note first that the BRLOs are the products of gas dynamics, rather 
than radiative transfer. RRL observers are conditioned to think about pressure broadening 
whenever they see broad line profiles or wide wings on the lines. We can rule out pressure 
broadening, however, as the source of the wide lines in BRLOs. The very strong dependence 
of the pressure broadened line width on principle quantum number, n (AVocn^) means that, 
if pressure broadening produces the wide lines, their widths should change dramatically (by 
factors of 5-400) over the range of n values observed for our BRLOs. In all cases for which 
there are observations of Hna transitions over a range of n, the observed changes in the 
width and/or strength of the broad recombination line component are much smaller than 
those predicted by pressure broadening. 

5.2.1. Champagne Flows 

Champagne models ( Tenor io-Tagle 1979, Bodenheimer, Tenorio-Tagle, & Yorke 1979) 
consider the evolution of young H II regions as they break out of the dense neutral cores 
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in which they form into a lower density diffuse medium. During the breakout phase, 
the ionized gas accelerates into the low density medium at speeds up to 30-50 km s^^. 
Subsequent analysis shows that this kind of acceleration also occurs in H II regions that 
form in clouds with power law density distributions, if the density gradients are sufficiently 
steep (Franco, Tenorio-Tagle, & Bodenheimer 1990). A champagne phase should be a 
common occurence among UC H II regions, if they really form as a result of ionization of 
the surrounding molecular core, because of their extremely large internal pressures. 

5.2.2. Ionization of Neutral Winds 

Hypersonic bipolar outflows of molecular gas are a common feature of protostellar 
sources, both in low mass star forming regions and in regions containing low and high 
mass pre main sequence objects (Bachiller 1996; Plume et al. 1992). These outflows have 
velocities ranging from a few up to about 100 km s^^. External ionization of molecular flows 
near O stars embedded in a young cluster could reproduce at least some of the features of 
BRLOs. Alternatively, if the turn-on time for the ionizing star is short enough and occurs 
soon enough after the outflow phase, the O star could ionize its own molecular flow from 
the inside. 

K3-50A is the only one of the broad line sources that is large enough (0.2x1.2 pc) to 
have a reasonably long crossing time (~2xl0^ years). In principle, then, this source would 
be a good candidate for the ionization of a fossil flow. In such a picture, we would expect 
the ionizing photons to move through the cavity evacuated by a fast molecular wind and 
strike the swept-up material in the cavity walls. The dense, UV-illuminated walls would 
appear as an elongated, limb-brightened shell. Existing continuum maps do not show such 
a shell (De Free et al. 1994). 
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5.2.3. Wind-Driven Bubbles 

The fast (1000-2000 km s^^) stellar winds from young O stars will interact strongly 
with the surrounding interstellar gas. After a short time (~10^ years, Shull 1980, 1982), the 
wind will create a bubble with a dense shell of swept-up material. If this bubble forms in a 
flattened structure, it may result in formation of a bipolar outflow of ionized material (Koo 
& McKee 1992). 

5.2.4^. Mass-Loaded Winds 

Redman, Williams, & Dyson (1995, 1996; Dyson, Williams & Redman 1995) and 
Lizano et al. (1996) have developed a scheme to produce long-lasting UC H II regions in 
which material photoevaporated from dense circumstellar clumps loads the O star wind and 
produces a dense outflow of ionized material bounded by a recombination front. 

5.2.5. Disk Winds 

Hollenbach et al. (1994) and Yorke and collaborators (Yorke & Welz 1993, 1996; 
Richling & Yorke 1997) have suggested that evaporative flows from massive circumstellar 
disks may offer a solution to the lifetime problem for UC H II regions. In their models, 
stellar UV photons ionize material in the disk and produce steady mass loss for periods 
of up to 10^ years. These flows are inherently bipolar. It is possible, however, to produce 
asymmetric regions by invoking external photoevaporation of the disks of nearby companion 
stars (Kessel, Yorke, & Richling 1998). Most of the papers discuss the terminal velocity of 
this evaporated material but do not explicitly model it. Kessel et al. (1998) present line 
proflles for optical transitions showing typical linewidths on the order of 40 km s^^. 



-24- 



5.2.6. Remarks on the Models 

Since the widespread existence of BRLOs was unknown at the time the various 
theoretical papers were pubhshed, the models do not explicitly address the issues raised 
by the nature of these sources. Nevertheless, the papers do allow us to categorize the 
potential problems of the various models. None of the models has an obvious explanation 
for the relatively narrow range in observed line widths. All but the last two models suffer 
from difficulties with the lifetime. They produce high velocity emission only for the order 
of a single crossing time and therefore cannot account for the number and frequency of 
observed BRLOs. We therefore consider only the mass loaded wind and photoevaporated 
disk models in detail. 

The mass-loaded wind models can only account for the bipolar morphology of most 
BRLOs if the underlying protostellar wind is strongly bipolar (the photoevaporated disk 
models, by contrast, only work if all of the BRLOs turn out to be bipolar). Mass loaded 
wind models also fail to explain the r~^ or steeper density gradients inferred from the 
continuum spectral indices. Mass loading should result in density profiles less steep than 
r~^. The existing versions of the mass-loaded wind and photoevaporated disk models also 
do not predict the very broad recombination line profiles we observe. 

On balance, the photoevaporated disk model has fewer problems than the mass-loaded 
wind model. There are observations that provide some supporting observational evidence 
for the presence of disks, if not for their involvement in the generation of the broad lines. 
For the best-studied BRLO, MWC 349, the hot dust surrounding the star shows a disk-like 
distribution with a size of <50xl00 AU (Leinert 1986). A similar size for the neutral disk 
has also been inferred from the high angular resolution radio continuum map at 1.3 cm 
(Mart in- Pint ado et al. 1993). Direct evidence of the interface region between the neutral 
disk and the ionized gas has been measured from the H30q; recombination line maser spikes 
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which also indicates a face-on disk of ~ 80 AU (Planesas, Mart in- Pint ado & Serabyn 
1992). The most detailed knowledge of the kinematics of the circumstellar gas comes from 
the IR recombination lines (Hamann & Simon 1986) and the observations of the HSOa 
maser line. Modeling of the IR line profiles indicates the presence of a neutral rotating 
disk and a bipolar ionized wind. The location and the radial velocities of the HSOa maser 
spikes are consistent with a keplerian disk rotating around a 30 Mq object (Planesas et 
al. 1992; Thum, Mart in- Pint ado & Bachiller 1992). The presence of a disk in S106 is also 
inferred from the morphology of the ionized gas (see Felli et al. 1985). Unfortunately, for 
the other BRLOs, the properties of the disks are unknown. We note, however, that these 
sources are excellent candidates for the detection of recombination line masers at millimeter 
wavelengths, similar to those observed in MWC 349 (Martin-Pintado et al. 1989). Detection 
of maser emission might allow us to observe the interface region between the neutral disk 
and the ionized material. 

5.3. Conclusions 

1) BRLOs are common in regions containing embedded OB stars. These objects have 
recombination line widths of 60-80 km s~^ and continuum spectra that rise as z/~°'^ to z/'"^'^ 
into the millimeter band. Most spatially resolved BRLOs have bipolar morphology. 

2) The short crossing times and high frequency of occurence of BRLOs implies both 
that these objects last for much longer than their expansion times and that many UGH II 
regions must begin their lives as BRLOs. 

3) The high continuum opacity of BRLOs at short centimeter wavelengths means that 
estimates of dust extinction based on near-IR lines and radio continuum measurements 
probably overestimate the extinction. 
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4) None of the published models of UCH II regions explicitly deals with all of the 
observed properties of BRLOs. If a photoevaporated disk is the source of the ionized bipolar 
emission region, the results imply that disks around massive young stars are both common 
and long-lived. 

This work was supported in part by National Science Foundation grant 95-30695 to 
the University of Texas at Austin. 
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Figure Captions 

Fig. 1. — Comparison of fits of one and two Gaussians to the H41a line toward 
NGC 7538 IRSl. (left) Single Gaussian fit to the region within ±15 km s"^ of the peak 
intensity and residuals to that fit. (right) Two Gaussian fit to the same profile (line 
parameters in Table ^. Most of the power in the residuals to this fit is due to either 
He41a or to the molecular feature seen at higher velocities (CH3CN). 

Fig. 2. — Comparison of fits of one and two Gaussians to the H41a line toward W49N. 
(upper left) Single Gaussian fit to the region within ±15 km s^^ of the peak intensity, (lower 
left) Residuals to the single Gaussian fit. (upper right) Two Gaussian fit to the same profile 
(see Table 0). (lower right) Residuals to the two Gaussian fit. The solid line shows a two 
Gaussian fit to the He41a line made using the widths of the hydrogen features and a constant 
velocity offset of -122.2 km s~^. 

Fig. 3. — Comparison of H39a (thin histogram) and H41a (dashed histogram) profiles in 
W49N. (upper panel) The dotted line shows the third feature in a three Gaussian fit to H39a;. 
For this fit, the width and position of the other two components were fixed to the values 
derived from a fit of two Gaussians to the H41q; profile, (lower panel) Overlay of the H39a 
profile and an H41q; profile to which the third component has been artificially added. 

Fig. 4. — Comparison of fits of one and two Gaussians to the H41a line toward G34.3±0.2. 
(upper left) Single Gaussian fit to the region within ±15 km s^^ of the peak intensity, (lower 
left) Residuals to the single Gaussian fit. (upper right) Two Gaussian fit. (lower right) 
Residuals to the two Gaussian fit and a two Gaussian fit to the He41a line (see text). 

Fig. 5. — Overlay of the H39a (solid line) and H41a (dashed line) spectra toward G34.3±0.2. 
The dashed line shows the sum of two Gaussians representing a compromise between fits 
to the two hydrogen transitions (V7v= 53.0 km s~^, AVAr=32.3 km s""^; Vb=41.6 km s~l, 



-33- 



-i> 



AVb=75.0 km s 

Fig. 6. — (iGft) H39a toward Orion H2 Peak 1. The smooth curve shows a single Gaussian 
fit to the observed spectrum (see Table 0). The residuals to this fit are below the profile. 
(right) H41a toward Orion IRc2. 

Fig. 7. — Recombination line spectra toward GlO.6-0.4. (top) H41a. The smooth curve 
shows a single Gaussian fit to the hydrogen line and a fit for the intensity of the He41a line 
with the width of this line fixed to the hydrogen line width and the velocity fixed to -122.2 
km s~^ with respect to the hydrogen line, {bottom) Fit to the H35a; line and He35a line. 

Fig. 8. — Single Gaussian fit to the H39a line toward DR21. The He39a line was fit holding 
the width equal to the H39q; width and the velocity at -122.2 km s~^ with respect to the 
hydrogen line. 
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Table 1. Observing Log 



Source 


R.A.(1950) 


Decl.(1950) 


Date 


Lines 


Orion H2 Peak 1 


05 32 46.3 


-05 23 54 


2/92 


H30a, H36a, H39a 


Orion IRc2 


05 32 47.0 


-05 24 24 


2/92 


H36a, H41a 


GlO.6-0.4 


18 07 30.6 


-19 56 28 


3/93 


H35a, H41a 


G34.3+0.2 


18 50 46.1 


1 11 10 


2/92 
3/93 


H36a, H39a, H41a 
H30a, H31a, H35a, H39a, H41a 


W49N 


19 07 49.8 


9 01 17 


2/92 
3/93 


H36a, H39a, H41a 
H30a, H31a, H35a, H41a 


S106 


20 25 32.8 


37 12 44 


2/92 
3/93 


H36a, H39a 
H35a, H41a 


DR21 


20 37 14.1 


42 08 55 


2/92 


H30a, H36a, H39a 




20 37 14.2 


42 09 06 


3/93 


H35a, H41a 


NGC 7538 


23 11 36.5 


61 11 49 


2/92 
3/93 


H30a, H36a, H39a 
H31a, H35a, H39a, H41a 
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Table 2. Continuum and recombination Line Parameters 



Source 


Line 




Tl 


V 


AV 


/TdV 


s.-) 


Te* ^''> 








K 


km s^^ 


km s^^ 


K km s-i 


Jy 


K 


Orion H2 Peak 1 


H39a 




0.27 


-4.34 


25.6 


7.3(0.25) 






Orion IRc2 


H41a 




0.54 


-3.57 


22.7 


13.0(0.3) 






GlO.6-0.4 


H41a 




0.67 


0.2(0.1) 


28.0(0.1) 


20.5(0.15) 


4.0 


9300 


G34.3+0.2 


H41a 


N 


0.45 


53.0 


32.3 


14.9(0.4) 


3.0 


9500 






B 


0.25 


41.6 


75.0 


21.1(0.6) 


4.2 


9500 


W49N 


H41a 


N 


0.77 


7.7 


30.6 


25.1 


6.0 


10500 






B 


0.18 


12.7 


70.8 


13.3 


3.0 


10500 


S106 


H41a 


N 


0.136 


-3.1(0.6) 


29.7(2.3) 


4.3(0.7) 


1.0 


10000 






B 


0.066 


-1.6(2.2) 


75(8) 


5.2(0.8) 


1.1 


10000 


DR21 


H41a 




1.24 


-0.84(0.06) 


32.6(0.1) 


43.2(0.2) 


10.6 


11000 


NGC 7538 


H41a 


N 


0.094 


-65.6(0.7) 


24.5(1.0) 


2.7(0.3) 


0.8 


12800 






B 


0.077 


-47.8(1.9) 


69.3(3.0) 


5.8(0.4) 


1.7 


12800 




H39a 


N 


0.095 


-65.6(0.6) 


26.8(1.9) 


2.7(0.3) 










B 


0.066 


-46.8(2.1) 


77.1(3.9) 


5.4(0.3) 







''The continuum flux densities for the normal (N) and broad (B) spectral components have 
been derived from the electron temperatures and the integrated line intensities 

'^LTE electron temperature derived from the ratio between the total integrated line intensity 
(B+N) and the measured continuum emission 



-36- 



Table 3. Summary of the properties of the broad recombination hne sources 



Source 


Morphology 


Spectral 


Size 


AV 


Cross, time 


References to 






index 


pc 


knis^^ 


102 yr 


col. 2 to 5 


Sgr B2F* 


Compact 


-1.5 


0.011 


59 


2 


a 


G5. 89-0.39 


Shell 


0.6, 2.2 


0.044 


64 


6 


b, c, d 


G34.3+0.2C 


Cometary 


0.9, -0.1 


0.035 


75 


5 


b, e 


W49N 


Unresolved 


0.6, 0.9 


0.07 


71 


10 


f, e 


K3-50A 


Bipolar 


0.5 


0.2x1.2 


66 


30-180 


g 


S106 


Bipolar 


0.73, -0.1 


0.07x0.2 


75 


9-12 


h,e 


MWC349 


Bipolar 


0.65 


0.001x0.03 


60 


0.2-0.6 


i,j,k 


NGC7538 IRSl 


Bipolar 


0.9 


0.003x0.02 


64 


0.5-3 


1, e 



* 4 sources, a) De Free et al. (1996), b) Wood & ChurchwcU (1989), c) See text, d) Zijlstra et al. 
(1990), e) this paper, f) De Free et al. (1996), g) De Free et al. (1994), h) Bally, Snell & Frcdmore 
(1983), i) Martin-Fintado et al. (1993), j) Harvey, Thronson & Gatley (1979), k) Hartmann, Jaffe & 
Huchra, (1980), 1) Gaume et al. 1995) 
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